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Bose-Einstein condensation in a two-dimensional, trapped, interacting gas 
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We study Bose-Einstein condensation phenomenon in a two-dimensional (2D) system of bosons 
subjected to an harmonic oscillator type confining potential. The interaction among the 2D bosons 
is described by a delta-function in configuration space. Solving the Gross-Pitaevskii equation within 
the two-fluid model we calculate the condensate fraction, ground state energy, and specific heat of 
the system. Our results indicate that interacting bosons have similar behavior to those of an ideal 
system for weak interactions. 
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The observation of the Bose-Einstein condensation 
(BEC) phenomenon in dilute atomic gases Jl|-|| has 
caused a lot of attention, because it provides opportu- 
nities to study the thermodynamics of weakly interact- 
ing systems in a controlled way. The condensate clouds 
obtained in the experiments consist of a finite number 
of atoms (ranging from several thousands to several mil- 
lions), and are confined in an externally applied confining 
potentials. The ground state properties of the condensed 
gases, including the finite size effects on the temperature 
dependence of the condensate fraction, are of primary 
interest. At zero temperature, the mean- field approxi- 
mation provided by the Gross-Pitaevskii equation [ 5| de- 
scribes the condensate rather well and at finite temper- 
atures a self-consistent Hartree-Fock-Bogoliubov (HFB) 
approximation is developed. || Path Integral Monte 
Carlo (PIMC) simulations on three-dimensional, in- 
teracting bosons appropriate to the current experimen- 
tal conditions demonstrate the effectiveness of the mean- 
field type approaches. Various aspects of the mean- field 
theory, as well as detailed calculations corresponding to 
the available experimental conditions are discussed by 
Giorgini et al. || 

In this work, we examine the possibility of BEC in 
a two-dimensional (2D) interacting atomic gas, under a 
trap potential. Such a system may be realized by making 
one dimension of the trap very narrow so that the oscil- 
lator states are largely separated. Possible experimental 
configurations in spin polarized hydrogen and magnetic 
waveguides are currently under discussion. Q The study 
of 2D systems is also interesting theoretically, since even 
though the homogeneous system of 2D bosons do not 
undergo BEC, Jl(| number of examples have indi- 
cated such a possibility upon the inclusion of confining 
potentials. We employ the two-fluid, mean-field model 



developed by Minguzzi et al. |12| to study the 2D Bose 
gas. Similar approaches |l3| are gaining attention be- 
cause of their simple and intuitive content which also pro- 
vide semi-analytical expressions for the density distribu- 
tion of the condensate. Recently, Mullin jl4|] considered 
the self-consistent mean-field theory of 2D Bose particles 
interacting via a contact interaction within the Popov 
and semi-classical approximations. His conclusions were 



that a phase transition occurs for a 2D Bose system, in 
the thermodynamic limit, at some critical temperature, 
but not necessarily to a Bose-Einstein condensed state. 
However, in the current experiments the finite number 
of atoms N prevent various divergences to give rise to 
behavior akin to non-interacting systems. 

Our work is motivated by the success of mean field, 
two- fluid models |12 1^] vis a vis more involved calcu- 
lations and direct comparison with experiments. In the 
following we briefly describe the two-fluid model of Min- 
guzzi et al. Jl^] and present our results for the 2D Bose 
gas. 

The condensate wave function \l/(r) is described by the 
Gross-Pitaevskii (GP) equation || 



2m 



V 2 *(r) + V cxt {r)^(r) + 2gnx{r)^{r) + <?* 3 (r) 

= M*(r) , (1) 



where g is the repulsive, short-range interaction 
strength, V r cx t(r) = mus 2 r 2 /2 is the confining (or trap) 
potential, n\ (r) is the distribution function for the non- 
condensed particles. We note that unlike in a three- 
dimensional system, g in our case is not simply related to 
the s— wave scattering length, but will be treated as a pa- 
rameter. In the two-fluid model developed by Minguzzi et 
al. H the non-condensed particles are treated as bosons 
in an effective potential V e g (r) = T4 x t(r) + 2gni(r) + 
2gi]/ 2 (r), and having the same chemical potential /x with 
that of the condensate. The density distribution is given 
by 



ni(r) 



1 



d 2 p 

(2nh) 2 exp{(p 2 /2m + V eS (r) - n)/k B T} - 1 ' 

(2) 



and the chemical potential is fixed by the relation 

p{E)dE 



N = N + 



exp [(E - n)/k B T] - 1 



(3) 



where N = J ^> 2 (r)d 2 r, is the number of condensed 
atoms, and the semi-classical density of states is calcu- 
lated using lOIPlPI 
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FIG. 1. The condensate fraction No/N as a function of 
T/T , for a system of N = 10 5 atoms. The various interaction 
strengths are described by the parameter r\ 



p(E) = 



d J r. 



(4) 



V Bti (r)<E 



The GP equation admits a simple solution within the 
Thomas-Fermi approximation, i.e. when the kinetic en- 
ergy term is neglected, 

* 2 (r) = -[»- V cxt (r) - 2gn 1 {r)]0(jji - V cxt (r) - 2< ? n 1 (r)) , 



fJ 



(5) 



where 9{x) is the unit step-function. Thomas-Fermi ap- 
proximation is regarded to be rather good except for the 
region close to the phase transition. |17| Minguzzi et al. 
|L2| have numerically solved the above set of equations 
self-consistently. They have also introduced a simpler ap- 
proximation scheme which treats the interaction effects 
perturbatively. Encouraged by the success of even the 
zero-order solution in describing the fully numerical self- 
consistent solution in the 3D case, we attempt to look at 
the situation in 2D. In a similar vein, we treat the interac- 
tions among the non-condensed particles perturbatively. 
To zero order in gni(r), the number of condensed parti- 
cles are calculated to be 



gm V Tilj / ' 
and the density of states is obtained as 

E/(huj) 2 if /i< , 

IHAE) { 2{E- li)/{Tiuj) 2 if -It, • /•. ■:// > 0) 
E/(hu) 2 if 2/i < E (n> 0) 



(6) 




FIG. 2. The temperature dependence of the chemical po- 
tential /x for various interaction strengths. 



If we use the above form of the density of states, valid 
for E > 0, then we obtain 



N = N a + t 2 



7T" 

T 



dilog ( 1 



,-a/t 



(8) 



where t = ksT /Tilo, and a = fj,/tko. The chemical poten- 
tial /i(iV, T) is obtained as the solution of this transcen- 
dental equation. 

In Fig. 1 we show the temperature dependence of the 
condensate fraction Nq/N for a system of N = 10 5 par- 
ticles, and for various values of the interaction strength. 
Also shown for comparison is the result for an ideal 
2D Bose gas in an harmonic trap, given by Nq/N = 
1 - (T/T ) 2 where k B T = hcj(N/({2))^ 2 . We observe 
that BEC-like behavior occurs for small values of the pa- 
rameter rj = mg/irh 2 , i.e. the weakly interacting system. 
Here we identify the BEC with the macroscopic occupa- 
tion of the ground state at T — and the depletion of it 
above Tq. As the strength of interactions is increased we 
find that the temperature dependence of N /N deviates 
from the non-interacting case more noticeably. Mullin 
Jl4j has argued that there is no BEC in 2D in the thermo- 
dynamic limit. We consider a system with finite number 
of particles, and we were able to obtain a self-consistent 
solution for the chemical potential for various values of 
the interaction strength as displayed in Fig. 2. We next 
evaluate the temperature dependence of the internal en- 
ergy (E) = [{E) nc (N - N )/2 + {E) c ]/N which consists 
of contributions from the non-condensed particles 



^o(f) 172 [^(2 C (3) 



(7) 



and the condensed particles 



-a/t 



(9) 
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FIG. 3. The ground state energy (i?) of the 2Dbosons as a 
function of temperature for various interaction strengths. The 
Maxwell-Boltzmann result is shown by the thin solid line. 



1/2 



a 



(10) 



In the above expressions C,( n ) is t ne Riemann zeta- 
function. The kinetic energy of the condensed particles 
is neglected in accordance with our Thomas-Fermi ap- 
proximation to the GP equation. In Fig. 3, we display 
the temperature dependence of (E) for different values 
of the interaction strength. The non-interacting energy 
is simply (E)/Nk B T = [C(3)/C(2)] (T/T ) 5 . For small 
rj, and T < To, the behavior of (E) resembles to that 
in a 3D system. As r\ increases, a bump in (E) devel- 
ops for T < To, which perhaps indicates the breakdown 
of the present approximation or an artifact of the cal- 
culation. We have no physical explanation for this be- 
havior. The corresponding results for the specific heat 
CV = d(E)/dT, are shown in Fig. 4. In contrast to the 
non-interacting case where a sharp peak at T = To is 
seen, the effects of short-range interactions smoothes out 
the transition. However, this smoothing is partly due to 
the finite number of particles in the system. || The ef- 
fects of interactions and finite number of particles are not 
disentangled in our treatment. 

It is a straightforward generalization to include the ef- 
fects of anisotropy within the present formalism. For 
an external potential of the type V cxt (r) = muj^(x 2 + 
A 2 y 2 )/2, where A = u> y /u) x is the anisotropy parameter, 
both Nq and po(E) depend inversely on A. Similarly, our 
analysis may be extended to study other power law po- 
tentials such as V ~ r 1 for which )«1 appears to be 
interesting. 
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Our calculations using the two- fluid model of Minguzzi 
et al. jL2| show that the BEC, in the sense of macro- 
scopic occupation of the ground state, may occur in a 2D 
trapped Bose gas when the short-range interparticle in- 
teractions are not too strong. As the interaction strength 
increases we could not find self-consistent solutions to the 



FIG. 4. The specific heat CV = d(E)/dT as a function of 
temperature for various interaction strengths. 



mean-field equations signaling the breakdown of our ap- 
proach. We note that instead of using the lowest order 
perturbation approach adopted here, the full solution to 
the self-consistent equations may alleviate the situation. 



Given the unclear nature of the phase transition [114 18 
in 2D and the interest of future experiments, we think 
it is worthwhile to perform first-principles calculations. 
Recent PIMC simulations |18| on 2D hard-core bosons 
confirm the possibility of BEC in the sense that a sharp 
drop in Nq/N around ksT c ~ 0.78 iV 1 / 2 is observed for 
finite systems. 

In summary, we have applied the mean-field, semi- 
classical two-fluid model for trapped interacting Bose 
gases to the case in two-dimensions. We have found that 
for a range of interaction strength parameters the behav- 
ior of the thermodynamic quantities resembles that of 
non-interacting bosons in a harmonic trap. 
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